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Abstract 

We examine classical Bogolyubov's model of a particle coupled to a heat bath which 
consists of infinitely many stochastic oscillators. Bogolyubov's result [I] suggests that, in 
the stochastic limit, the model exhibits convergence to thermodynamical equilibrium. It has 
recently been shown that the system does attain the equilibrium if the coupling constant is 
small enough [T5]. We show that in the case of the large coupling constant the distribution 
function ps(q,p,t) — > pointwise as t ^ 00. This implies that if there is convergence to 
equilibrium, then the limit measure has no finite momenta. Besides, the probability to find 
the particle in any finite domain of phase space tends to zero. This is also true for domains 
in the coordinate space and in the momentum space. 



1 Introduction 

If two bodies with different temperatures are in contact, they will eventually have the same tem- 
perature. The inverse process of "temperature separation" does not occur if we do not act on the 
system by anything. This phenomenon is referred to as irreversibility. It seems paradoxical since 
equations of mechanics (Newton's equation) and quantum mechanics (Schrocdinger's equation) are 
time-reversible. This problem has been discussed for a long time and a lot of outstanding scien- 
tists such as Boltzmann, Poincare, Gibbs, Birkhoff, Bogolyubov and others tried to solve it. As a 
result, new approaches and techniques have been developed PQ-[B]. One of the recently developed 
techniques is a stochastic limit (see [5] and references therein). 
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The idea of Bogolyubov's model pQ, that considers the behavior of one particular oscillator 
under the action of many other stochastic oscillators, was later further developed [7] . The quantum 
analogue of Bogolyubov's model has been studied in details as well (see [H]-[II] and references 
therein). 

In this paper we first briefly describe Bogolyubov's model. Bogolyubov pQ suggested a toy 
model that could represent a system in contact with a thermostat. The thermostat is modelled by 
an infinite number of oscillators whose initial coordinates and momenta are random variables with 
thermal (Gibbs) distribution. The system is represented by a single oscillator whose coordinate 
and momentum are arbitrarily fixed at the initial instant. The system interacts with the thermo- 
stat with some coupling constant. It is expected that asymptotically the system gets the same 
temperature as the thermostat, i.e. the coordinate and momentum of the single oscillator will 
obey the Gibbs distribution. In his paper [T] Bogolyubov proved an estimate of the distribution 
function ps(q,P, t) in some interval of t, which suggests that, in the stochastic limit [5], the model 
exhibits convergence to thermodynamical equilibrium. Bogolyubov's model is simple enough to 
prove theorems or make explicit calculations in some particular cases. It has recently been shown 
that the system does attain the equilibrium if the coupling constant is small enough [12) . 

In this paper we examine the model in the case of the large coupling constant. We find that in 
this case the distribution function ps(q,p,t) — > pointwise as t — > oo. This implies that if there 
is convergence to equilibrium, then the limit measure has no finite momenta and is not the Gibbs 
function. Besides, the probability to find the particle in any finite domain of phase space tends to 
zero. This is also true for domains in the coordinate space and in the momentum space. 

The outline of the paper is as follows. In Sec. 2 we formulate a mathematical model and set 
out Bogolyubov's results. In Sec. 3 we give a theorem about attaining equilibrium in a particular 
case of small coupling constants [12]. And in Sec. 4 we consider another particular case of large 
coupling constants. In Sec. 5 we discuss the results. 



2 Model and Bogolyubov's results 

The Hamiltonian and Hamilton equations. The following model is considered. There is 
an oscillator (the system) and a set of N oscillators (the thermostat) with the following total 
Hamiltonian: 

1 1 N N 

n— 1 n— 1 

where p, q, u> and p n , q n , uj n are momenta, coordinates and frequencies of the first oscillator 
and those of the set of oscillators, respectively; e and a n are positive numbers and play a role of 
coupling constants. In what follows we imply e talking about a small or large coupling constant. 
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The corresponding Hamilton equations are 

^ 2 ^n . 2 dq n / n \ ,-, / n \ 

-^T + u n q„ = -ea n q, p n = —, p n (0) = P n , q n (0) = Q n , 

d 2 q , 2 ^ dq (2) 



dt 2 , uj 2 q = -e^2a n q n , P = ~^' P(0) = Po, 9(0) = qo- 

n=l 



The model parameters a n , w n , P n , Q n , po, qo satisfy the following conditions. The initial 
momentum and coordinate of the system p , q are arbitrary real numbers: p , qo € M. 

The parameters a n and the frequencies u n satisfy the conditions corresponding to transition 
to a continuous spectrum as N — > oo: 

2 ru 2 roo 

E ^ - / J ^ E ^ - / J M d - (3) 

for V i/ > 0. J{v) is a continuous positive function and / J{y)dv < oo. 

Jo 

The initial momenta and coordinates of the set of oscillators (the thermostat) P n and Q n are 
random variables with the distribution function 

p(( n ,9 n ) = exp ^ - ^ f; + W y n )j (4) 



such that 



/ p(Cn,fln)dCl---dCjV<Wl--.dfliV = l, 



where f €l and fc, T are positive numbers. Physically, fc and T are Boltzmann constant and 
temperature, respectively. 

Bogolubov's results. Let us introduce new variables E n and (p n as follows: 

Qn = -cos<^„, P n = -\ / 2E n sin<p„, (5) 

so that _E„ = ^(-Pn + ^nOn) arc initial energies. Further, let 

jv . 

^w = E^ !E ^. ( 6 ) 

n=i w " 



iV 



y/*2P 

/Ar(t) = -E a ™ - cos (w„* + <£n) (7) 



n=l W " 
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and Ujv(i) be a solution of the integro-differential equation 

(8) 



v'jf(t) + uj 2 v N (t) = e 2 I K N (t - t)v(t)cIt, 



v N (0)=0, v' N (0) = l. 
Then the solution q(t),p(t) of equations @ reads [T] 



q(t) = q v' N (t) + p v N (t) +e v N (t- T)f N {r)dT, 



o, 



p{t) = goUjv(*) +Po^(*) +£ v' N (t- T)f N (r)d 



(9) 



Note that the dependance of the solutions q(t) and on A~ is implied. 

Bogolyubov pQ showed that as N — > oo the solution Ujv(t) along with its first and second 
derivatives converges uniformly in any finite interval to v(t). The latter is a solution of the 
following integro-differential equation: 



v"(t) + uj 2 v(t) = e 2 f Q(t - T)v'{ T )dT, 
Jo 



(10) 



«(0)=0, t/(0) = l, 
where 

Q(t)= / J(v)(l-cosvi)dv. 
Jo 

According to Bogolubov [1] we can formulate 

Theorem 1. There exists a limit of the probability density of random values q(t), p(t) for any 
t > as N — > oo: 

p S (t,q,p) = *(q - q*(t),p- p*(t),t). 
The limit is meant in the following sense: 

pa2 r b2 

m ProMa-i < a(t) < aoM < v(t) < M = I 

TV- 
Here 

q*(t) = q v'(t) + PQ v(t), p*{t) = q Q v"(t) +p v'{t) (12) 

and 

=. r + 1 ( C^ 2 -2B^ + A V 2 \ 

The coefficients A = A(t), B = B(t) and C = C(t) are derived from the identity 



pa2 pb2 

lim Prob{ ai < q(t) < a 2 ,h < p(t) < b 2 } = / p s (t,Z,rj)d£dr). 

N ^°° Jal Jbl 



(11) 



A(t)\ 2 + 2B(t)Xfi + C{t)ii 2 = e 2 kT / J(v) 



{Xv(x) + ^v'{x)}e~ lvx dx 



dv. (14) 
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From the identity (|2"5)) it is clear that, first, A > and, second, AC — B 2 > 0. The latter is obvious, 
because the right hand side is positive for any A and /i, hence, the discriminant 4(B 2 — AC) < 0. 

The second important result in PQ yields an estimate of the limit function ps(t,q,p) in some 
interval with respect to t and is formulated as 

Theorem 2. For Ve > 0, V/3 > a > 0, and for any sequence {At e } such that At £ — > oo, 

e 2 At e — > as e — > we have for Vt € [ 



t+At e 



7TT~ / 



(15) 



where cr(e) — > as e — > 0, and is some explicit expression which tends to the Gibbs distribution 
with temperature T as t — > oo. 

However, Theorem 2 tells us nothing at all about the asymptotic behavior of ps(t,q,p) as 
t — > oo. In [T2] a particular case has been considered and some kind of an asymptotics, which 
tends to the Gibbs function, has been found. The respective theorem is formulated in the next 
section. 



3 Particular case with a small coupling constant e 

We shall leave clarifying what should be considered as a small or a large coupling constant till the 
next section. 

Theorem 3. Let J{v) 6 C(R) n Li(M) be an even rational function, and all its critical points 
in C are of the first order. Then for any a > there is that for any e: < |e| < £o there exists 
such to (e) that when t > to (e) we have for any p, q € M 

lo ( E + Eoe- 25 ^ - 2^fEE ~e- s ^ t cos ((lu + e 2 Imp)t + fo - ip)\ 
Ps{q ' P > t) - 2nkT(l-e- 2 ^) - eXP { (1 - e -^)fcT ' < ^ 



(16) 



where 8(e) and p(e) are determined by the function J(v). Besides, 

V2E 



cosip, p = — y/2E$ixi(p > 

LO 

q = cost^o, po = —V2E sax(po, 

LO 

p z LO z q z 
i.e. E = — H — is the energy. 

From Theorem 3 one can easily see that the asymptotics as t — > oo tends to the Gibbs function. 
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4 Particular case with a large coupling constant e 



Let us consider a particular case with 

J{v)-- 



1 



a>0, 6>0. 



(17) 



a + bv 2 

Obviously, this function satisfies the Theorem 3 conditions. 

Proposition 1. If J(v) is taken as specified in (jTTJ) and function v(t) in integro-differential 



equation (|10| is triply continuously differentiable then (|10j) takes the form of a third-order differ- 
ential equation: 



v'"(t) + Q 1/2 v"(t)+^v'(t)+^Q 
v{0) = 0, v'(0) = 1, u"(0) = 0, < t < +c» 



1/2 „ £ 2 7T 

a; 

26 



v(t) = 0, 



Proof. First of all, let us calculate Q(t) with residues. 

1 — cos vt 



Q(t)= / J(v) (1 - cos vt)dv = 



1 



, , dv = — 
a + bv 2 2 



1 - e l 



bv 2 



dv = 



2biyfa/b 2bi^/a/b 



2Vab 



l-e 



Integrating the right-hand side of (|10l) by parts we can write it in the following form: 



v"{t) + uj 2 v{t) = e 2 / v{r)Q'{t - r)d 



Taking into account the explicit formula we obtain: 



'(t)+Jv{t) = — F(t), 



2b 



where 



F(t) 



v(t) exp 



-(t-r) )dT. 



Clearly, F(t) satisfies the equation: 



dF 
~dt 



F + v{t). 



(18) 



(19) 



(20) 



(21) 



(22) 



Then we differentiate both parts of equation (|2T|) and obtain the third-order differential equation 
in (fT5|) . One more initial condition, which is the value of the second derivative v"(0), directly 
comes from (fTU|) if we let t = 0. The proposition is proved. 



G 



The corresponding characteristic equation for (1181) is 



X* + J^-X 2 +lj 2 X+ J-uj z --- = 0, (23) 




or 



(24) 



At this point we can formulate the difference between a small and a large coupling constant. 
If equation (|23j) has two complex roots, which differ by order of e 2 from iu> and — iui (the roots 
are purely imaginary in the case of £ = 0), and one real root which differs by the same or- 
der from — y/a/b, then this is the case of a small coupling constant. And this is the case of a 
large coupling constant when (|23|) has three real roots: two negative and one positive. We can 
make sure that the characteristic equation can have two negative and one positive roots. Let 
uo 2 = l/3,£7r/26 = 4 and a/b = 9. Then the characteristic equation (|2"3"|) takes the form: 

X 2 + l) (A + 3) =4. 



3, 

It is easy to check that the last equation has three real roots whose approximate values are 
— Ai w —2.2723, — A2 ~ —1.5691 and A3 0.8414. In the case of a large coupling constant we shall 
prove the following 

Theorem 4. Let equation (|2"3l have three real roots two of which are negative and one is 
positive: — Ai, — A 2 and A3, where Ai > 0, A2 > 0, A3 > 0. Further, let A3 < A 2 , A3 < Ai and 
A3 < yj a/b. Then the pointwise limit with respect to q and p is equal to zero: 

lim p s (t,q, P ) = 0. (25) 

t — >+oo 

Proof. In order to prove the theorem we explicitly calculate A, B and C. According to the 
conditions of the theorem the solution of equation ([TH]) is 

v(t) = C ie - Xlt + C 2 e- X2t + C 3 e X3t . (26) 
From the initial conditions we have: 

d = 



A3 — A2 


(A 2 


-Ai)(AH 


-A 3 )' 




Ai — A3 




(A 2 


-Ai)(A.H 


-A 3 )' 




Ai + A 2 





°2 = T, TTTT^v (27) 



C 3 = 



(A 2 + A 3 )(Ai + A 3 )' 
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Then we find A(t), B(t) and C(t) from the equality: 

roc rt 

A{t)\ 2 + 2B(t)Xft + C(t)n 2 = s 2 kT / J{v) / {Xv(x) + iiv'{x)}e- ivx dx 

Jo Jo 

where J(u) 



dv, (28) 



1 



a + bv 2 ' 

Let us introduce and Si. i= 1,6, as follows: 



/ {Xv(x) + nv'(x)}e- lux dx 
Jo 



h+h + h + h + h + h. 



(29) 



/>OG 

/ J{") 
Jo 



where 



{Xv(x) + fiv'(x)}t 



Si 



dv = S 1 + S 2 + S 3 + S 4 + S 5 + S 6 , (30) 
J(y)lidv. 



Straightforward, but tedious calculations give (some intermediate calculations are carried out 
in Appendix A) 

; (A-MAi) 2 



h = c{- 



X\ + v 2 



(l-2e" Alt cos^ + e" 2Alt ) 



-2A 2 i\ 



J 2 = Cl { \ 2 /iA2 2 } (1 - 2e- x * cosvt + , 
X 2 + v 

(l-2e A3t cos^ + e 2A3t ). 



T _^ 2 (A + M 3 ) 2 

A3 + V 1 



I4 — 2ClC2 



(A - fj,Xi)(\ - llX 2 ) 



(l + e -(Ai+A a )t)( AlA2 + 1/ 2)_ 



( A 2 + ^ )(A 2 + I/2) 

-(e~ Al * + e" A2t )(AiA 2 + v 2 )cosvt - ^(e~ Alt - e~ A2 *)(Ai - A 2 )sinz/t] 



(31) 
(32) 
(33) 

(34) 



T _ or r (A + //A 3 )(A- M A 2 ) 
/5 - 2C2C3 (A| + ^)(A| + ^) 



J/ 2 ) l [ 



= (A3-A 2 )t Vl/ 2 



)K-A 3 A 2 )- 



(35) 



-(e Aat + e~ A2t )(^ 2 - A 3 A 2 )cosz/t + j/(e A3 * - e~ A2 *)(A 2 + X 3 )smut] 



r _ r, n n (X - nXi)(X + (iX 3 ) 
(Af + z/2)(A 2 + z/ 2 ) 



( 1 + e (A3-A 1 )t )(l/ 2_ AiAa) _ 



(36) 



-{er Xlt + e Xst )(v 2 - AiA 3 )cosj/i- ^(e~ Alt - e A3 *)(A 3 + Ai)sini/i] . 
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Si 



s 2 = - 2 



7rC?(A- mAi) 2 
2(a-6A 2 ) 



^C 2 (A-/zA 2 ) 2 



(l + e- 2Alt ) 



2(a-6A 2 ) 
_ ^7 2 (A + M A 3 ) 2 



(1 + e" 2 * 2 *) 



Ai 



A 2 



2e 



2e 



-Ait 



-A 2 t 



-Ait 



Ai 

,-A 2 t 

"aT" 

-A 3 t 



3 2(a - 6A§) 

_ 7tC 1 C 2 (A-mAi)(A-mA 2 ) 



(a-6A 2 )(a-6A 2 ) 

6 3/2 



^ _ e -(Ai+A 2 )t^ 



i+a 2 )A 2a-6(A 2 + A 2 ) 



Ai + A 2 



+ 



+ A 2 )t _ -(Xl + y/f)t _ -(A 2 + ^/f )t 



- be'^ (e- A2t - e- Alt ) (A 2 - Ai) 



S 5 = 



ttC 2 C 3 (X - AiA 2 )(A 4- //A 3 ) 



(a-6A 2 )(a-6A 2 ) 



, X3-A 2 )t>! 2a - KAj + A 2 ) 
A 2 — A3 



&3/2 

A 2 A 3 ^ + Va&) (1 + e 



(A 3 -A 2 )t _ -(A 2 +-/f)t _ (Ag-vT)* 



6e"V** (e- A2t - e A3t ) (A 2 + A3) 



^C 1 C 3 (A- A tA 1 )(A + A tA 3 ) 
(a-6A 2 )(a-6A 2 ) 

& 3/2 



ri_ P (As-Ai)A 2a-&(A 2 + A 2 ) 
V ) Ai — A 3 



_ (AiAs^ + V^) (l + e (A3 - Al) * - e-^+x/B* - e^"^') 

- be-Vt* ( e ~ Xlt - e A3t ) (Ai + A3) 
We define Pi, i = 1, 6 as follows: 

5! = (A- M Ai) 2 Pi, 

52 = (A — M 2 ) 2 P 2 , 

5 3 = (A + A tA 3 ) 2 P 3 , 

5 4 = (A - juAi)(A - ^A 2 )P 4 , 

5 5 = (A- M A 2 )(A + mA 3 )P 5 , 

5 6 = (A- A tA 1 )(A + A tA 3 )P 6 . 



(37) 

(38) 
(39) 



(40) 



(41) 



(42) 



(43) 
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Now we use the theorem conditions A3 < A2, A3 < Ai, A3 < v/ a/b. 

Then all Si but S3 tend to constants. The latter grows exponentially: S3 oc e 2A3t . (Hereafter, 
the notation cx has the same meaning as in S3 = const ■ e 2Xat + o(e 2A3t ) as t — > 00). 

From definition we have: 
A 



e 2 kT 



= P 1 +P 2 + P 3 + P 4 + P 5 + P 6 , 



L A 1 P 1 -A 2 P 2 + A3P3-^(Ai + A 2 )F4-kA 2 -A 3 )P5-^(Ai-A3)P 6 , (44) 



C 



e 2 kT 



A 2 Pi + A 2 P 2 + A 2 P 3 + AiA 2 P 4 - A 2 A 3 P 5 - AiA 3 P 6 , 



(45) 



where P3 oc e 2A;,t and the other Pi oc const as t — > +00. 

Then, we investigate the behavior of AC — B 2 . It is clear that the terms quadratic in P3 are 
eliminated. The question is whether the coefficient in front of terms linear in P3 is zero or not. 

Making necessary substitutions from (|44p we obtain: 

= - A 2 ) 2 P 4 2 - J(A 2 + A 3 ) 2 P 5 2 - i(A! + A 3 ) 2 P 2 + 

+(Ai - A 2 ) 2 PP 2 + (Ax + A 3 ) 2 PiP 3 + (Ai + A 3 ) (A! - A 2 )PiP 5 + 
+ (A 2 + A 3 ) 2 P 2 P 3 + (A 2 + A 3 )(A 2 - A!)P 2 P 6 + (A 2 + A 3 )(Ax + A 3 )P 3 P 4 + 

+^(Ai - A 2 )(A 2 + A 3 )P 4 P 5 + i(A 2 - Ai)(Ai + A 3 )P 4 P 6 - ^(Xi + A 3 )(A 2 + A 3 )P 5 P 6 . 

From the latter expression we obtain that the behavior is as follows as t — ► +00: 

AC- B 2 oc aP 3 , 
a = (Ai + A 3 ) 2 Pi + (A 2 + A 3 ) 2 P 2 + (A 2 + A 3 )(Ai + A 3 )P 4 . 

Proposition 2. lim a > 0. 

t — >oc 

Proof. The behaviors of P 4 , P 2 , P 4 are: 

irCf ( 1 [b\ 



(46) 



P oc 



2(a-6A 2 ) I Ai 



«"^(j5i;)' (47) 

(a - 6Af)(a - b\%) \ Ai + A 2 
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The Viet theorem for characteristic equation (|23j) takes the form: 



A3 — Ai — A2 — — 



A1A2 — A2A3 — A1A3 — U) 
\ \ \ - £2?r f°- 2 

AlA 2 A 3 — — a -U! . 

2b V 



(48) 



An intermediate result is lim a = (3R (the Viet theorem (|4"5|) is already partially used and 

t — >oo 

relations (|27p are taken into account), where 

'(A3-A2) 2 , (A3-AO 2 4(A 3 -A a )(A 3 -A 1 ) 



Ai 



A 2 



Ai + A2 



A l/i \ \2 , A i/\ \ \2 „(Af + Af)(A 3 - A 2 )(A 3 - Ai) 
— (A 3 -A 2 J +— (A3-A1) -2 — — 

Ai A 2 Ai + A2 



2(a-6A 2 )(a-6A2)(A 2 -Ai) 2 
Using Viet relations (|4*5|) again we have: 

(A 2 - Ai) 2 



R = a 



AiA 2 (Ai + A 2 ) 



AiA 2 + -J + 



-b 



A 2 + Aj (Az-AQ 2 
2 AiA 2 (Ai + A 2 ) 



AiA 2 + T ) +b 



" \ , (A 2 -A 1 ) 2 (A 2 + A 1 ) 2 



2AiA 2 (A! + A 2 ) 



Hence, 



lim a = — — — — — > 0, 

t^oo 6AiA 2 (Ai + A 2 ) 



which proves the proposition. 

Calculating the exponent in (1131) . The exponent we calculate is: 

Cf 2 - 25£r? + Ai] 2 



where 
and 



11 2(AC~B 2 ) ' 

£,=q-q*{t), r)=p-p*(t) 



q*(t) = q v'(t) +p v(t), p*(t) = q Q v"(t) +p v'(t). 

As t — > co, g*(t) and p*(i) behave as follows: 

g^txCge^teoAg+Po) 
p*(t)oc A 3 C3e A3t ( 9 oA 3 +po)- 



(49) 



(50) 



(51) 



(52) 



(53) 



(54) 
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Hence, for arbitrary finite p and q 



It yields 



7?(i) oc -A 3 C3e A3t (goA3 +Po)- 



n = C^-2Bfr ? + V k (goA3+P o) 2 C3 2 e^^(c - 2i?A 3 + AXj) 



2{AC-B 2 ) 
Using (gl]) and(g6|) we obtain: 



2(AC-B 2 ) 



2aXi ( 1 

lim n = n - 31 

t— » + oo 



•A, ' = const. 



And finally, 



lim ps(t,q,p)\p,q= c <mst = lim - q*(t),p - p*(t), t) = lim 



t — >+oo 



-A 3 t 



i^ +00 2nyjAC - B 2 



2a\ 2 ( 1 



which proves the theorem. 

Corollary 1. Letp p (p,i)= / p s (q,p,t)dq and p q {q,t) = / p s (q,p,t)dp. Then 
lim p p (p,t) = lim p q (q,t) = 0. 

i — >oo i — >oo 

Proof. From the explicit expression (fT3|) one can easily calculate 



P P (p,t) 



1 



and 



Then equations (J3H) and ([55)1 yield: 



1 



exp 



V2^A 



exp 



lim p p (p,t) = lim = 



exp 



(p-p*W) 2 

2C 

(g-g*ffl) 2 

2A 



A|C|(goA 3 +Po)e^ 

^3 



since P3 ex e 2A3t as t — > 00. Analogously, lim p q (q,i) = 0. 
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5 Discussion 



Let us calculate the mean coordinate (q) and momentum (p), and their standard deviations. Using 
Corollary 1 we obtain 



(q) = / O Pq (O,t)d0 = q*(t), 



( P ) = / d Pp (6,t)de = p*(t), 



((q-q*(t)) 2 ) = (6-q*(t)) 2 Pq (6,t)d6 = A(t), ((p-p*(t)) 2 ) = (6-p*(t)) 2 p p (6,t)dO = C(t) 

(62) 

Again, taking relations (|44[) and ([55)) into account we see that the behavior of the mean values 
and the standard deviations is exponential: 

(?) oc C 3 e A3t ( go A3 +po), (p) oc XsCse^iqoXs+Po), 

(63) 

{(q-q*(t)) 2 ) cxe 2A3t , {(p-p*(t)) 2 ) oce 2A3 *. 

This behavior seems strange. Since qo and po are arbitrary real numbers, (q) and (p) may tend 
either to the positive or negative infinity depending on sign{qo\j, + po). It appears that the particle 
goes away to the infinity exponentially, and its standard deviation increases exponentially as well. 
However, this strange behavior is explained by 

Theorem 5. Characteristic equation (|2"3"|) has a positive root if, and only if, the Hamiltonian |T|) 
is not positive-definite as a quadratic form of (2N + 1) variables (q, qx, . . . , qw,Pi, ■ ■ ■ ,Pn) as 
iV — ► oo. 

Proof. We use the Silvester criterion to find out when the quadratic form H(q, q%, . . . ,qN,Pi, ■ ■ ■ ■, 
from ([T]) is positive-definite. Its doubled matrix is 

/ 



L> 2 


ea\ 


ECti 


£Q!jv 





. 


■ M 


EOL\ 


CO 2 











. 


. 


ea 2 





u% 








. 


. 














. 


. 


ECIN 








... u% 





. 


. 














1 


. 


. 














1 . 


. 
















. 

















. 


■ 1 / 



(64) 



V 

In this case it is sufficient to consider only first (iV + 1) determinants in the Silvester criterion. 
The n-th determinant D n can be easily calculated and is as follows: 



D„ 



col 



' ^ io 2 

1=1 1 



(65) 
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From the latter formula one can see that it is sufficient to require only Dn to be positive: then 
the rest determinants are positive. If Dn > then 

.2 W 2 



^ < -R—- (66) 

,,,2 



= 1 



Hence, as N — > oo (fMjl turns into 



uj 2 2\fbaoj 2 ,„„. 
£ 2 < — = — • (67) 

J{r)dT n 



Then the right-hand side of ([241 (which is equivalent to ([25]) ) is less than u 2 J — . Then from the 



plot of the left-hand side of (as a function of A) it is clear that if (f6"T)) is true, equation 
cannot have a positive root. And vice versa, if (f67|) is not satisfied eq. (|23|) has a positive root, 
but the Hamiltonian is not positive-definite. Similar divergencies associated with non-positivity of 
the density matrix were found in the quantum analogue of Bogolyubov's model [11] . 

It is worth noting that the exponential runaway of the particle mean coordinate and momentum 
is not intrinsic to the stochastic character of the thermal bath oscillators. In the deterministic case 
(when one solves © with certain initial data) this also may occur. Indeed, in the simplest case 
when E n — (or, equivalently, P n — and Q n — 0) it is easy to notice from © (in this case 
/jv(i) = 0) that lim q(t) = q*(t), and as we have already seen q*{t) cx C^X^qa + po)e Xst . 

N— »oo 



6 Conclusion 

It is possible for any coupling constant e find such a, b and u> that the limit (the limit N — > oo is 
computed) distribution function tends to zero as t — > +oo and p and q are fixed. This implies that 
if there is convergence to equilibrium, then the limit measure has no finite momenta. Moreover, 
the probability to find the particle in any finite domain of the phase, coordinate or momentum 
space tends to zero, although the integral all over the space equals to 1. This phenomenon might 
be related to the fact that, as it follows from Theorem 5, the Hamiltonian is not positive-definite 
in this regime for large N. 
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A Intermediate calculations of Si, i = 4, 5, 6 

S 4 = CiC 2 (A - /iAi)(A - mA 2 ) [(1 + e-( Al+A2)t )S 4 ,i - (e~ Alt + e"^ 4 )^- 



where 



-(e- A2t -e-^*)(A 2 -A 1 )5 4 , 3 ] , 



AiA 2 + v 2 



/oo ^ \ 

_ JM (ATT^)(Ai + , 2 ) 

+ v 2 ) cos vt 



dv, 



a f°° il \ (^1^2 - 



S^ 3 = / J{y) 



5*4,1 — 



{a - bX\){a - bX 2 ) 



2 )(Ai+^) 

sin i4 

(A? +I /2)(Ai + l /2) 

\2 , \2 



2a-b(X 2 + X 2 ) y— 
v \ — - Vab + AiA; 



Ai + A 2 



63/2 



'4.2 



(a-6A?)(o-6A|) 



a(e- Alt + e-^ 1 ) - b(Af e~ Aat + Aj e - Alt ) 
Ai + A 2 



+ 



53/2 

+ I AiA 2 ^ - Vab) e'Vft 



'4.3 



a(e- Alt - e- A2 *) 6(Aie- Alt - A? e" A2 *) " 



(a - bX 2 )(a - bX 2 ) 

7rCiC 2 (A- A1 A 1 )(A- A iA 2 ) 



A 2 - Af 



A 2 _ A l 



S 4 = 



(a-6A?)(a-&A|) 

6 3/2 



1+A2)t x 2a-b(X 2 + X 2 ) 



^ _ e -(A 1+ A 2 )^ ^-0 



+ A 2 



+ 



+ ^A 2 ^ - V^) (l + e -(Ai+A 2 )t _ e -(A 1+> /f)* _ e -(A 2 +V¥)^ _ 
- fce-VT* (e- A2t - e - A i*) (A 2 - Ai)" . 



(68) 

(69) 
(70) 
(71) 
(72) 

(73) 
(74) 



(75) 
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S 5 = C 2 C 3 (X - M A 2 )(A + M A 3 ) [(1 + e (A3 - A2) *)^ 5 ,i - (^ A2t + e Xat )S 5 , 2 - 



(e -A 2 t_ e A 3 t )(A2 + A3)5 I 



where 



— A2A3 



roc v -2 . 
^5,1 = / J ( V )7V2~, : 

J -co (A| + v 



2 )(A| + ^ 2 ) 



5*5,2 — 



'5,3 



/OO 
"CO 



J(iz) 



(i^ 2 — A2A3) COS ft 

(Af + ^XAi+^2) 

f sin i4 
(AI + ^ 2 )(A| + ^) 



dv, 



dv. 



5*5,1 — 



55,2 — 



(a-b\l){a-b\l) 



7T 



(a-&A 2 )(a-6A 2 ) 

- ( Vab + A2A3 



6(A 2 + A 3 ) - ( \/a6 + A 2 A 3 ^ 

v a 



e -A 3 t_ e -A 2 t A 2 e -A 2 t _ A 2 e -A 3 t 

+ 6- 



A3 — A 2 
53/2 



'5,3 



(a-b\l){a-b\l) 



A3 — A 2 



A3 - A 2 . 



A3 - A 2 . 



7rC 2 C3(A-^A 2 )(A + /xA 3 ) 



(a-6A 2 )(a-6A 2 ) 



(1 _ P (A 3 -A 2 )A 2a-&(A 2 + A 2 ) 
V J A 2 -A 3 



A 2 A 3 - 



^ + V^j (l + e (*s-*0* - e -(^+>/f)* - e (^-^)t) 



- be - ^* (e- A2t - e A3t ) (A 2 + A 3 ) 

One can see that 56 is obtained from S5 by substitution 2 — > 1: 
7rC 1 C 3 (A-Mi)(A + MA3) 



5 6 = 



(a-6A 2 )(a-6A 2 ) 

& 3/2 



/ (As-aoA 2a-&(A 2 + A 2 ) 
V 1 J A! -A 3 



(A 3 -Ai)t _ e -(Ai+Vf )* _ e ( A3 -V f ' 



- (AiA 3 ^ + \/^) (l + e 



0" 



(76) 

(77) 
(78) 
(79) 
(80) 

(81) 
(82) 



(83) 



(84) 
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